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CLAIMS 



[Claim(s)] 

[Claim 1] To the subject-copy image data which it comes to arrange in the shape of a 
matrix, the pixel value of the 1st number It is an identification information embedding 
method for embedding the identification information which consists of a numerical signal 
of the 2nd number below said 1st number. Relate with said each numerical signal the 
combination of the basis function which intersects perpendicularly mutually, and the 2nd 
carries out number generation. For every combination of said basis function by calculating 
total of the product of the value and the brightness value of a pixel of each basis function to 
the location of each pixel in subject-copy image data which intersects perpendicularly The 
weighting factor of the 2nd number corresponding to the combination of the basis function 
of the 2nd number is computed, respectively. The embedded function which is the number 
of multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take while including the value which said weighting factor can take in a domain 
for said every numerical signal is referred to. The output value among two or more input 
values which are in agreement with the numeric value of the numerical signal [ from ] 
Specify the input value nearest to the weighting factor corresponding to the combination of 
the basis function related with the numerical signal, and all the weighting factors of said 
2nd number so that it may become said input value specified about each, and the same 
value The identification information embedding method to the image data characterized by 
changing the pixel value of said image data. 

[Claim 2] From the processing-object image data where said identification information was 
embedded by the identification information embedding method to image data according to 
claim 1 Are the identification information extract approach for extracting said 
identification information, and relate with said each numerical signal the combination of 
the basis function which intersects perpendicularly mutually, and the 2nd carries out 
number generation. For every combination of said basis function by calculating total of the 
product of the value and the brightness value of a pixel of each basis function to the 
location of each pixel in processing-object image data which intersects perpendicularly The 
weighting factor of the 2nd number corresponding to the combination of the basis function 
of the 2nd number is computed, respectively. The embedded function which is the number 
of multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take while including the value which said weighting factor can take in a domain 
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for said every weighting factor is referred to. The identification information extract 
approach from image data that the identification information characterized by computing 
said embedded value of a function to the weighting factor was embedded. 
[Claim 3] To the subject-copy image data which it comes to arrange in the shape of a 
matrix, the pixel value of the 1st number Are an identification information embedding 
method for embedding the identification information which consists of a numerical signal 
of the 2nd number below said 1st number, and orthogonal transformation is performed to 
each pixel value of said subject-copy image data. The multiplier distribution data which 
come to arrange the weighting factor of the 1st number in the shape of a matrix are 
generated. The weighting factor of the 2nd number chosen from each weighting factor 
which constitutes this multiplier distribution data Make said which numerical signal 
correspond and the embedded function which is the number of multi-pair Ichinoseki which 
includes in range the numeric value which said numerical signal can take while including 
the value which said weighting factor can take in a domain is referred to for each [ these ] 
numerical signal of every, respectively. The input value nearest to the weighting factor in 
said multiplier distribution data corresponding to the numerical signal is specified from 
from among two or more input values whose output value of the corresponds with the 
numeric value of the numerical signal. The identification information embedding method 
to the image data characterized by performing reverse orthogonal transformation to the 
multiplier distribution data with which the weighting factor concerned in said multiplier 
distribution data was permuted, and the permutation of the weighting factor 
corresponding to all numerical signals was made with the specified input value. 
[Claim 4] From the processing- object image data where said identification information was 
embedded by the identification information embedding method to image data according to 
claim 3 Are the identification information extract approach for extracting said 
identification information, and orthogonal transformation is performed to each pixel value 
of said processing-object image data. The multiplier distribution data which come to 
arrange the weighting factor of the 1st number in the shape of a matrix are generated. The 
weighting factor of the 2nd number corresponding to said each numerical signal is taken 
out from this multiplier distribution data. The embedded function which is the number of 
multi-pair Ichinoseki which includes the numeric value which said numerical signal can 
take while including the taken-out value which said weighting factor can take for every 
weighting factor in a domain in range is referred to. The identification information extract 
approach from image data that the identification information characterized by computing 
said embedded value of a function to the weighting factor was embedded. 
[Claim 5] To the subject-copy image data which it comes to arrange in the shape of a 
matrix, the pixel value of the 1st number Are an identification information embedding 
method for embedding the identification information which consists of a numerical signal 
of the 2nd number below said 1st number, and a two-dimensional discrete cosine transform 
is given to each pixel value of said subject-copy image data. The frequency distribution 
data which come to arrange the weighting factor of the 1st number in the shape of a matrix 
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are generated. The weighting factor of the 2nd number chosen from each weighting factor 
which constitutes this frequency distribution data Make said which numerical signal 
correspond and the embedded function which is the number of multi-pair Ichinoseki which 
includes in range the numeric value which said numerical signal can take while including 
the value which said weighting factor can take in a domain is referred to for each [ these ] 
numerical signal of every, respectively. The input value nearest to the weighting factor in 
said frequency distribution data corresponding to the numerical signal is specified from 
from among two or more input values whose output value of the corresponds with the 
numeric value of the numerical signal. The identification information embedding method 
to the image data characterized by giving a two-dimensional reverse discrete cosine 
transform to the frequency distribution data with which the weighting factor concerned in 
said frequency distribution data was permuted, and the permutation of the weighting 
factor corresponding to all numerical signals was made with the specified input value. 
[Claim 6] From the processing-object image data where said identification information was 
embedded by the identification information embedding method to image data according to 
claim 5 Are the identification information extract approach for extracting said 
identification information, and a two-dimensional discrete cosine transform is given to 
each pixel value of said processing-object image data. The frequency distribution data 
which come to arrange the weighting factor of the 1st number in the shape of a matrix are 
generated. The weighting factor of the 2nd number corresponding to said each numerical 
signal is taken out from this frequency distribution data. The embedded function which is 
the number of multi-pair Ichinoseki which includes the numeric value which said 
numerical signal can take while including the taken- out value which said weighting factor 
can take for every weighting factor in a domain in range is referred to. The identification 
information extract approach from image data that the identification information 
characterized by computing said embedded value of a function to the weighting factor was 
embedded. 

[Claim 7] Said embedded function is an identification information embedding method to 
image data given in any of claims 1, 3, and 5 characterized by being a periodic function 
they are. 

[Claim 8] Said embedded function is the identification information extract approach from 
image data that identification information given in any of claims 2, 4, and 6 characterized 
by being a periodic function they are was embedded. 

[Claim 9] Said embedded function is an identification information embedding method to 
image data given in any of claims 1, 3, and 5 characterized by being the continuous 
periodic function they are. 

[Claim 10] Said embedded function is the identification information extract approach from 
image data that identification information given in any of claims 2, 4, and 6 characterized 
by being the continuous periodic function they are was embedded. 

[Claim 11] Said embedded function is an identification information embedding method to 
image data given in any of claims 1, 3, and 5 characterized by it being narrow when 
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spacing of two or more input values which take the same output value has a small input 
value, and being large when an input value is large they are. 

[Claim 12] Said embedded function is the identification information extract approach from 
image data that identification information given in any of claims 2, 4, and 6 characterized 
by it being narrow when spacing of two or more input values which take the same output 
value has a small input value, and being large when an input value is large they are was 
embedded. 

[Claim 13] To the subject-copy image data which it comes to arrange in the shape of a 
matrix, the pixel value of the 1st number It is identification information embedded 
equipment for embedding the identification information which consists of a numerical 
signal of the 2nd number below said 1st number. An embedded function maintenance 
means to hold the embedded function which is the number of multi-pair Ichinoseki which 
includes in range the numeric value which said numerical signal can take while including 
the value which said weighting factor can take in a domain, An orthogonal transformation 
means to generate the multiplier distribution data which perform orthogonal 
transformation to each pixel value of said subject-copy image data, and come to arrange 
the weighting factor of the 1st number in the shape of a matrix, While making said which 
numerical signal correspond, the weighting factor of the 2nd number chosen from each 
weighting factor which constitutes this multiplier distribution data, respectively The input 
value nearest to the weighting factor in said multiplier distribution data corresponding to 
the numerical signal is specified from from for each [ these ] numerical signal of every 
among two or more input values whose output values of said embedded function 
correspond with the numeric value of the numerical signal. A weighting-factor 
permutation means by which the specified input value permutes the weighting factor 
concerned in said multiplier distribution data, Identification information embedded 
equipment to the image data characterized by having a reverse orthogonal transformation 
means to perform reverse orthogonal transformation to the multiplier distribution data 
with which the permutation of a weighting factor was made with this weighting- factor 
permutation means. 

[Claim 14] From the processing-object image data where said identification information 
was embedded by the identification information embedded equipment to image data 
according to claim 13 An embedded function maintenance means to hold the embedded 
function which is an identification information extractor for extracting said identification 
information, and is the number of multi-pair Ichinoseki which includes in range the 
numeric value which said numerical signal can take while including the value which said 
weighting factor can take in a domain, An orthogonal transformation means to generate 
the multiplier distribution data which perform orthogonal transformation to each pixel 
value of said processing-object image data, and come to arrange the weighting factor of the 
1st number in the shape of a matrix, The identification information extractor from image 
data with which the identification information characterized by having the fetch means 
which takes out the weighting factor of the 2nd number corresponding to said each 
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numerical signal from this multiplier distribution data, and a taken out calculation means 
to compute said embedded value of a function for every weighting factor was embedded. 
[Claim 15] Orthogonal transformation is made to perform to a computer to each pixel value 
of the subject-copy image data which come to arrange the pixel value of the 1st number in 
the shape of a matrix. The multiplier distribution data which come to arrange the 
weighting factor of the 1st number in the shape of a matrix are made to generate. The 
weighting factor of the 2nd number chosen from each weighting factor which constitutes 
this multiplier distribution data Which numerical signal of the identification information 
which consists of the 2nd numerical signal is made to correspond, respectively. The 
embedded function which is the number of multi-pair Ichinoseki which includes in range 
the numeric value which said numerical signal can take while including the value which 
said weighting factor can take in a domain for each [ these ] numerical signal of every is 
made to refer to. The input value nearest to the weighting factor in said multiplier 
distribution data corresponding to the numerical signal is made to specify from from 
among two or more input values whose output value of the corresponds with the numeric 
value of the numerical signal. The computer-readable medium which stored the program to 
which reverse orthogonal transformation is made to perform to the multiplier distribution 
data with which the weighting factor concerned in said multiplier distribution data was 
made to permute, and the permutation of the weighting factor corresponding to all the 
numerical signals of said identification information was made with the specified input 
value. 

[Claim 16] Orthogonal transformation is made to perform to a computer to each pixel value 
of the processing-object image data where the identification information which consists of a 
numerical signal of the 2nd number while coming to arrange the pixel value of the 1st 
number in the shape of a matrix is embedded. The multiplier distribution data which come 
to arrange the weighting factor of the 1st number in the shape of a matrix are made to 
generate. The weighting factor of the 2nd number corresponding to said each numerical 
signal is made to take out from this multiplier distribution data. The computer- readable 
medium which stored the program which makes the embedded function which is the 
number of multi-pair Ichinoseki which includes in range the numeric value which said 
numerical signal can take refer to while including the taken-out value which said 
weighting factor can take for every weighting factor in a domain, and makes said 
embedded value of a function to the weighting factor compute. 



DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[The technical field to which invention belongs] This invention relates to the technique 
which embeds the identification information which shows rightful claimants, such as a 
copyright person, into the digital image data which circulates through various media, such 
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as a communication line, and the technique of extracting identification information from 

such digital image data. 

[0002] 

[Description of the Prior Art] With development of digital technique in recent years and 
development of multimedia society, various information is changed into digital data and 
circulates widely through various kinds of media, such as a communication network, 
satellite communication, and CD-ROM. The problem of the duplicate of digital data has 
followed on circulation of the digital data in such multimedia society unescapable. 
[0003] If the duplicate of this digital data is the lawful range, it can contribute to 
development of the culture in multimedia society, but since digital data can reproduce 
without degradation when it corresponds to an unauthorized use like [ in case it is used on 
direct commerce ], there is a possibility that profits loss of rightful claimants (an author, a 
copyright person, a publishing-right person, rightful claimant of neighboring right, etc.) 
may become serious. 

[0004] Therefore, the technique for considering as the proof at the time of this image data 
being unjustly reproduced about image data by embedding identification information into 
data in a mode with difficult recognition on appearance conventionally is proposed. For 
example, the following techniques are proposed in J.Cox et al."Secure Spread Spectrum 
Watermarking for Multimedia", NEC Reserch Institute, and Technical Report 95" 10. That 
is, according to this technique, orthogonal transformation of the image data is carried out, 
two or more weighting factors which correspond to each dot location of identification 
information among the weighting factors of each basis function obtained by conversion are 
chosen, the value of each dot of identification information is added to the value of each 
selected weighting factor, and reverse orthogonal transformation is performed to all the 
weighting factors containing the multiplier after addition, consequently identification 
information embedded settled image data is generated. 
[0005] 

[Problem(s) to be Solved by the Invention] However, since the weighting factor obtained by 
carrying out orthogonal transformation of the original image data (before identification 
information embedding) took various values according to the above-mentioned Prior art, it 
was impossible to have specified that to which the value of each dot of identification 
information is added from from among each weighting factor obtained by carrying out 
orthogonal transformation of the identification information embedded settled image data, 
therefore, when a duplicate object appears while saving and managing the original image 
data and identification information embedded settled image data, respectively in order to 
carry out the technique mentioned above At the same time it extracts the value of each dot 
of identification information by subtracting each weighting factor obtained from each 
weighting factor obtained by carrying out orthogonal transformation of the identification 
information embedded settled image saved and managed by carrying out orthogonal 
transformation of the original image data By subtracting each weighting factor obtained 
from each weighting factor obtained by carrying out orthogonal transformation of the 
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duplicate object by carrying out orthogonal transformation of the original image data, the 
value of each dot of identification information had to be extracted and the identity of both 
identification information had to be proved. 

[0006] Thus, since identification information embedded settled image data and 
subject-copy image data had to be saved and managed at the duplex according to the Prior 
art, while the certification at the time of data control or detection was complicated, as 
compared with the circulating amount of data, twice as many storage as this was needed. 
In the case of the database treating the image data of about 1000 sheets, or the database 
which must update data frequently like a newspaper database, such a problem is serious 
especially. 

[0007] In addition, in order to avoid the problem of duplex management of this data, 
permuting a part of weighting factor obtained by carrying out orthogonal transformation of 
the subject-copy image by the value of identification information itself is also considered, 
but such replacement has a possibility that identification information may be immediately 
recognized by the duplicate person only by carrying out orthogonal transformation of the 
image data, and rewriting of data may be made while worsening remarkably the image 
quality of the image after reverse orthogonal transformation. 

[0008] Then, the technical problem of this invention is the mode in which the identification 
information extract from identification information embedded settled image data is 
possible without subject-copy image data and which cannot be recognized to a duplicate 
person in view of the above problem. In the identification information extract approach 
from image data that identification information was embedded by the identification 
information embedding method and the equipment, this approach, or equipment to the 
image data which can embed identification information at image data, without degrading 
image image quality and equipment, and a list It is related with the computer-readable 
medium which stored the program as which a computer is operated as such identification 
information embedded equipment or an identification information extractor. 
[0009] 

[Means for Solving the Problem] Invention given [ each ] in a claim is made in order to 
solve the above-mentioned technical problem. The pixel value of the 1st number to namely, 
the subject-copy image data which it comes to arrange in the shape of a matrix as 
invention according to claim 1 was shown in the principle Fig. of drawing 1 It is an 
identification information embedding method for embedding the identification information 
which consists of a numerical signal of the 2nd number below said 1st number. Relate with 
said each numerical signal the combination of the basis function which intersects 
perpendicularly mutually, and the 2nd carries out number generation (Si). For every 
combination of said basis function by calculating total of the product of the value and the 
brightness value of a pixel of each basis function to the location of each pixel in 
subject-copy image data which intersects perpendicularly The weighting factor of the 2nd 
number corresponding to the combination of the basis function of the 2nd number is 
computed, respectively (S2). The embedded function which is the number of multi-pair 
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Ichinoseki which includes in range the numeric value which said numerical signal can take 
while including the value which said weighting factor can take in a domain for said every 
numerical signal is referred to. The output value among two or more input values which 
are in agreement with the numeric value of the numerical signal [ from ] Specify the input 
value nearest to the weighting factor corresponding to the combination of the basis 
function related with the numerical signal (S3), and all the weighting factors of said 2nd 
number so that it may become said input value specified about each, and the same value 
(S4) which changes the pixel value of said image data - it is characterized by things. 
[00 10] Thus, according to this invention, the numerical information of identification 
information is not embedded at a direct weighting factor, but the thing nearest to the value 
of an original weighting factor among two or more input values of the embedded function 
which makes this numerical information an output value is permuted by this weighting 
factor. Therefore, the difference of the value of the weighting factor before and behind a 
permutation is very small, and ends. Therefore, image quality of the image data after 
modification of a pixel value is not not much degraded. Moreover, since it corresponds with 
which numeric value of the numerical information which constitutes identification 
information with an embedded function, even if the weighting factor after modification of a 
pixel value does not have subject-copy image data, it can extract identification information 
from the image data after modification of a pixel value. Furthermore, if a third person does 
not know an embedded function even if he calculates a weighting factor, he cannot detect 
what kind of numeric value each weighting factor supports. Therefore, since a third person 
cannot get to know the contents of identification information, he cannot change this 
identification information. 

[0011] The pixel value of subject-copy image data may be the brightness value of 
monochrome image, the brightness value of each chrominance signal of an RGB code, or a 
brightness value of a YCC signal, and may be a color difference value of a YCC signal. 
[0012] The number of dots of subject -copy image data may be the same number in every 
direction, and may differ in all directions. The numeric value which each numerical signal 
of identification information can take may be binary, and may be a value beyond it. 
[0013] Each numerical signal of identification information may be image information 
which may be arranged together with the shape of Rhine and has been arranged in the 
shape of a matrix. Calculation of a weighting factor may be performed only to the 
combination of the basis function corresponding to each numerical signal of identification 
information, and it may be carried out to the combination of all basis functions according to 
orthogonal transformation. According to the former, the processing total of calculation 
processing of a weighting factor decreases. According to the latter, after permuting by the 
input value which the basis function corresponding to each numerical signal of 
identification information constructed [ input value ], and had the weighting factor which 
******** specified, the pixel value of image data can be changed only by carrying out 
reverse orthogonal transformation to all weighting factors, as this orthogonal 
transformation - two-dimensional - alienation - cosine conversion, two-dimensional 
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dispersion sign conversion, and a two-dimensional Hadamard transform can be used. 
[0014] The embedded function may be held with the gestalt of a function expression, and 
may be held with the gestalt of the table which specified the relation of one output value 
over two or more input values. Moreover, as long as it is the number of multi-pair 
Ichinoseki, an embedded function may be a periodic function or there may be at a periodic 
function. [ no ] If it is a periodic function, since a function expression will become simple, 
the equipment which carries out this invention can be simplified. That is, since there is a 
solution for every point that only the integral multiple of a period shifted those 
intersections if only it asks for an intersection with the numeric value of each numerical 
signal of the identification information in 1 period of an embedded function, specification of 
an input value becomes easy. In addition, even if a component with a bigger weighting 
factor changes the value of that weighting factor, the error by this modification stops easily 
being able to be visible when it thinks from a viewpoint of the image quality of the image 
data after modification of a pixel value subjectively. Moreover, the resistance over the 
alteration of various image data becomes strong, so that the degree of modification is large. 
What is necessary is here, just to enlarge the period of an embedded function, in order to 
enlarge the degree of modification. An embedded function with which a period becomes 
large may be used as the value of a weighting factor becomes large so that the value of a 
weighting factor may be changed more greatly and the value of a weighting factor may be 
smaller changed from this to a smaller weighting factor to a bigger weighting factor. 
[0015] Invention according to claim 2 from the processing-object image data where said 
identification information was embedded by the identification information embedding 
method to image data according to claim 1 Are the identification information extract 
approach for extracting said identification information, and relate with said each 
numerical signal the combination of the basis function which intersects perpendicularly 
mutually, and the 2nd carries out number generation. For every combination of said basis 
function by calculating total of the product of the value and the brightness value of a pixel 
of each basis function to the location of each pixel in processing-object image data which 
intersects perpendicularly The weighting factor of the 2nd number corresponding to the 
combination of the basis function of the 2nd number is computed, respectively. While 
including the value which said weighting factor can take in a domain for said every 
weighting factor, with reference to the embedded function which is the number of 
multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take, it is characterized by computing said embedded value of a function to the 
weighting factor. 

[0016] Invention according to claim 3 the pixel value of the 1st number to the subject-copy 
image data which it comes to arrange in the shape of a matrix Are an identification 
information embedding method for embedding the identification information which 
consists of a numerical signal of the 2nd number below said 1st number, and orthogonal 
transformation is performed to each pixel value of said subject-copy image data. The 
multiplier distribution data which come to arrange the weighting factor of the 1st number 
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in the shape of a matrix are generated. The weighting factor of the 2nd number chosen 
from each weighting factor which constitutes this multiplier distribution data Make said 
which numerical signal correspond and the embedded function which is the number of 
multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take while including the value which said weighting factor can take in a domain 
is referred to for each [ these ] numerical signal of every, respectively. The input value 
nearest to the weighting factor in said multiplier distribution data corresponding to the 
numerical signal is specified from from among two or more input values whose output 
value of the corresponds with the numeric value of the numerical signal. It is characterized 
by performing reverse orthogonal transformation to the multiplier distribution data with 
which the weighting factor concerned in said multiplier distribution data was permuted, 
and the permutation of the weighting factor corresponding to all numerical signals was 
made with the specified input value. 

[0017] Invention according to claim 4 from the processing-object image data where said 
identification information was embedded by the identification information embedding 
method to image data according to claim 3 Are the identification information extract 
approach for extracting said identification information, and orthogonal transformation is 
performed to each pixel value of said processing-object image data. The multiplier 
distribution data which come to arrange the weighting factor of the 1st number in the 
shape of a matrix are generated. The weighting factor of the 2nd number corresponding to 
said each numerical signal is taken out from this multiplier distribution data. While 
including the taken-out value which said weighting factor can take for every weighting 
factor in a domain, with reference to the embedded function which is the number of 
multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take, it is characterized by computing said embedded value of a function to the 
weighting factor. 

[0018] Invention according to claim 5 the pixel value of the 1st number to the subject-copy 
image data which it comes to arrange in the shape of a matrix Are an identification 
information embedding method for embedding the identification information which 
consists of a numerical signal of the 2nd number below said 1st number, and a 
two-dimensional discrete cosine transform is given to each pixel value of said subject-copy 
image data. The frequency distribution data which come to arrange the weighting factor of 
the 1st number in the shape of a matrix are generated. The weighting factor of the 2nd 
number chosen from each weighting factor which constitutes this frequency distribution 
data Make said which numerical signal correspond and the embedded function which is 
the number of multi-pair Ichinoseki which includes in range the numeric value which said 
numerical signal can take while including the value which said weighting factor can take 
in a domain is referred to for each [ these ] numerical signal of every, respectively. The 
input value nearest to the weighting factor in said frequency distribution data 
corresponding to the numerical signal is specified from from among two or more input 
values whose output value of the corresponds with the numeric value of the numerical 
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signal. It is characterized by giving a two-dimensional reverse discrete cosine transform to 
the frequency distribution data with which the weighting factor concerned in said 
frequency distribution data was permuted, and the permutation of the weighting factor 
corresponding to all numerical signals was made with the specified input value. 
[0019] Invention according to claim 6 from the processing-object image data where said 
identification information was embedded by the identification information embedding 
method to image data according to claim 5 Are the identification information extract 
approach for extracting said identification information, and a two-dimensional discrete 
cosine transform is given to each pixel value of said processing-object image data. The 
frequency distribution data which come to arrange the weighting factor of the 1st number 
in the shape of a matrix are generated. The weighting factor of the 2nd number 
corresponding to said each numerical signal is taken out from this frequency distribution 
data. While including the taken-out value which said weighting factor can take for every 
weighting factor in a domain, with reference to the embedded function which is the 
number of multi-pair Ichinoseki which includes in range the numeric value which said 
numerical signal can take, it is characterized by computing said embedded value of a 
function to the weighting factor. 

[0020] The embedded function which can be set they to be [ any of claims 1, 3, and 5 ] is a 
periodic function, and invention according to claim 7 is specified. The embedded function 
which can be set they to be [ any of claims 2, 4, and 6 ] is a periodic function, and invention 
according to claim 8 is specified. 

[0021] Invention according to claim 9 is that the embedded function which can be set they 
to be [ any of claims 1, 3, and 5 ] is a continuous periodic function, and is specified. 
Invention according to claim 10 is that the embedded function which can be set they to be 
[ any of claims 2, 4, and 6 ] is a continuous periodic function, and is specified. 
[0022] Invention according to claim 11 is specified because it is narrow when spacing of two 
or more input values which take the same output value of the embedded function which 
can be set they to be [ any of claims 1, 3, and 5 ] has a small input value, and it is large 
when an input value is large. 

[0023] Invention according to claim 12 is specified because it is narrow when spacing of two 
or more input values which take the same output value of the embedded function which 
can be set they to be [ any of claims 2, 4, and 6 ] has a small input value, and it is large 
when an input value is large. 

[0024] Invention according to claim 13 the pixel value of the 1st number to the subject-copy 
image data which it comes to arrange in the shape of a matrix It is identification 
information embedded equipment for embedding the identification information which 
consists of a numerical signal of the 2nd number below said 1st number. An embedded 
function maintenance means to hold the embedded function which is the number of 
multi-pair Ichinoseki which includes in range the numeric value which said numerical 
signal can take while including the value which said weighting factor can take in a domain, 
An orthogonal transformation means to generate the multiplier distribution data which 
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perform orthogonal transformation to each pixel value of said subject-copy image data, and 
come to arrange the weighting factor of the 1st number in the shape of a matrix, While 
making said which numerical signal correspond, the weighting factor of the 2nd number 
chosen from each weighting factor which constitutes this multiplier distribution data, 
respectively The input value nearest to the weighting factor in said multiplier distribution 
data corresponding to the numerical signal is specified from from for each [ these ] 
numerical signal of every among two or more input values whose output values of said 
embedded function correspond with the numeric value of the numerical signal. It is 
characterized by having a weighting-factor permutation means by which the specified 
input value permutes the weighting factor concerned in said multiplier distribution data, 
and a reverse orthogonal transformation means to perform reverse orthogonal 
transformation to the multiplier distribution data with which the permutation of a 
weighting factor was made with this weighting-factor permutation means. 
[0025] Invention according to claim 14 from the processing-object image data where said 
identification information was embedded by the identification information embedded 
equipment to image data according to claim 13 An embedded function maintenance means 
to hold the embedded function which is an identification information extractor for 
extracting said identification information, and is the number of multi-pair Ichinoseki 
which includes in range the numeric value which said numerical signal can take while 
including the value which said weighting factor can take in a domain, An orthogonal 
transformation means to generate the multiplier distribution data which perform 
orthogonal transformation to each pixel value of said processing-object image data, and 
come to arrange the weighting factor of the 1st number in the shape of a matrix, It is 
characterized by having the fetch means which takes out the weighting factor of the 2nd 
number corresponding to said each numerical signal from this multiplier distribution data, 
and a taken-out calculation means to compute said embedded value of a function for every 
weighting factor. 

[0026] Invention according to claim 15 makes orthogonal transformation perform to a 
computer to each pixel value of the subject-copy image data which come to arrange the 
pixel value of the 1st number in the shape of a matrix. The multiplier distribution data 
which come to arrange the weighting factor of the 1st number in the shape of a matrix are 
made to generate. The weighting factor of the 2nd number chosen from each weighting 
factor which constitutes this multiplier distribution data Which numerical signal of the 
identification information which consists of the 2nd numerical signal is made to correspond, 
respectively. The embedded function which is the number of multi-pair Ichinoseki which 
includes in range the numeric value which said numerical signal can take while including 
the value which said weighting factor can take in a domain for each [ these ] numerical 
signal of every is made to refer to. The input value nearest to the weighting factor in said 
multiplier distribution data corresponding to the numerical signal is made to specify from 
from among two or more input values whose output value of the corresponds with the 
numeric value of the numerical signal. It is characterized by being the computer-readable 
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medium which stored the program to which reverse orthogonal transformation is made to 
perform to the multiplier distribution data with which the weighting factor concerned in 
said multiplier distribution data was made to permute, and the permutation of the 
weighting factor corresponding to all the numerical signals of said identification 
information was made with the specified input value. 

[0027] Invention according to claim 16 makes orthogonal transformation perform to a 
computer to each pixel value of identification information ******** rare ****** 
processing-object image data which consists of a numerical signal of the 2nd number while 
coming to arrange the pixel value of the 1st number in the shape of a matrix. The 
multiplier distribution data which come to arrange the weighting factor of the 1st number 
in the shape of a matrix are made to generate. The weighting factor of the 2nd number 
corresponding to said each numerical signal is made to take out from this multiplier 
distribution data. The embedded function which is the number of multi-pair Ichinoseki 
which includes the numeric value which said numerical signal can take while including the 
taken-out value which said weighting factor can take for every weighting factor in a 
domain in range is made to refer to. It is characterized by being the computer-readable 
medium which stored the program which makes said embedded value of a function to the 
weighting factor compute. 
[0028] 

[Embodiment of the Invention] Hereafter, the gestalt of operation of this invention is 
explained based on a drawing, embedded [ which is an identification information 
embedding method to the image data based on this invention, and the gestalt of operation 
of equipment ] - business - the computer for an extract which is the gestalt of operation of 
the identification information extract approach from image data that a computer and 
identification information were embedded, and equipment It is constituted so that it may 
explain below, and all well-known kinds of orthogonal transformation processing and 
reverse orthogonal transformation processing can be used conventionally and the 
subject-copy image and identification information (however, identification information <= 
subject-copy image) of all magnitude can be processed. 

(Outline of an identification information embedded **** extract) First, before explaining 
the concrete configuration of this operation gestalt, the principle of the identification 
information embedding to the image data based on this operation gestalt and the outline of 
the identification information extract from image data in which identification information 
was embedded are explained. 

[0029] Here, the subject-copy image data for identification information embedded 
(monochrome image data or brightness data extracted from NTSC color picture data) 
presuppose that it consists of pixels of a NxN individual. And the identification 
information for embedded presupposes that it is data which consist of a signal of L 
(L<=NxM) individual with which the predetermined value was given, respectively. 
[0030] embedded - business - two or more blocks (however, the number of pixels of each 
block L or more pieces) which divided the whole subject-copy image data mentioned above 
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or subject- copy image data, and were acquired boil a computer, respectively, it receives, 
and performs orthogonal transformation. A subject-copy image is divided and orthogonal 
transformation is performed for extracting a signature image without degradation to each 
block, even when the image with which embedded processing was performed is cut off 
partially and reproduced. However, if magnitude of each partial image is made small too 
much, since the image with which embedded processing was performed will deteriorate, 
what is considered as 8x8 pixels or more is desirable. 

[0031] On the occasion of orthogonal transformation, L combination of the basis function 
which was associated to each signal of identification information and which intersects 
perpendicularly mutually is prepared. And for every combination of each basis function, 
total of the product of the value and brightness value of each basis function to the location 
of each pixel of subject-copy image data which intersects perpendicularly is calculated, and 
it asks for the weighting factor of the combination of each basis function. In addition, it 
says that the direction of the variable of each basis function lies at right angles mutually in 
the subject-copy image which consists of a pixel of the NxM individual "which intersects 
perpendicularly", or a block here. 

[0032] next, embedded - business - a computer embeds the value of each signal which 
constitutes identification information at the weighting factor of the combination of the 
basis function associated beforehand, respectively, this time - embedded - business - a 
computer does not add the value of each signal to a direct weighting factor, does not 
overwrite a weighting factor, but prepares a predetermined embedded function for every 
signal, and transposes the value corresponding to the value of each signal to the original 
weighting factor with this embedded function. 

[0033] The embedded function prepared for every signal may be a function of a definition 
which is different for every signal, respectively, and may be a function of a definition 
common to each signal. Even if it is which case, each embedded function includes the value 
which a weighting factor can take in the domain, and includes the value which each signal 
of identification information can take in the range. Moreover, each embedded function is 
the number of multi-pair Ichinoseki corresponding to an output value with two or more 
same input values included in a domain. Since processing will specifically be simplified if it 
is a simple function as this embedded function, a periodic function (when spacing of two or 
more input values which take the same output value has a small input value, it is narrow, 
and it is the continuous periodic function large when an input value is large) as shown, for 
example in the following formula (l) is desirable. 
[0034] 

[Equation l] 

[0035] embedded - business - a computer counts backward all the input values of the 
embedded function which constitutes identification information and which makes the 
value of the signal an output value for every signal. And the input value which specified 
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that to which a difference with the weighting factor of the combination of the basis function 
beforehand matched with the signal becomes the smallest, and was specified from from 
among all the computed input values is transposed to the value of the original weighting 
factor. 

[0036] after that - embedded - business a computer rewrites subject-copy image data or 
the value of each pixel of a block so that the value of the weighting factor of the 
combination of each basis function after orthogonal transformation may replace and it may 
become a next value. Thus, identification information embedded settled image data is 
obtained. In addition, since the variation by replacement of the value of a weighting factor 
is stopped to the minimum as mentioned above, degradation of identification information 
embedded settled image data is stopped by the minimum. Moreover, since the value of 
each signal of identification information cannot be known if an embedded function is not 
known even if a third person does orthogonal transformation of the identification 
information embedded settled image data, it is impossible to change this identification 
information. 

[0037] The computer for an extract performs orthogonal transformation to each of two or 
more blocks (however, the number of pixels of each block L or more pieces) which divided 
the whole processing-object image data or processing-object image data, and were acquired. 
At this time, L combination of the basis function which was associated to each signal of 
identification information and which intersects perpendicularly mutually is prepared like 
****** of the identification information to subject-copy image data. And for every 
combination of each basis function, total of the product of the value and brightness value of 
each basis function to the location of each pixel of processing- object image data which 
intersects perpendicularly is calculated, and the weighting factor of the combination of 
each basis function is called for. In addition, as these basis functions, the thing of ****** G f 
the identification information to subject-copy image data and the same thing are usually 
used. 

[0038] The computer for an extract inputs each called-for weighting factor into a 
corresponding embedded function, and calculates an output value. And each output value 
of the embedded function called for by doing in this way is arranged in as the order of a list 
of the signal of the identification information beforehand associated to the corresponding 
weighting factor. Then, if processing-object image data is identification information 
embedded settled image data, the output value put in order is in agreement with 
identification information. Thus, even if the computer for an extract does not have 
subject-copy image data, if only it holds each embedded function, it can extract 
identification information from identification information embedded settled image data. 
And this embedded function may be used in common to various subject-copy image data. 
Therefore, the total amount of the data which should be saved and managed decreases 
sharply as compared with the conventional thing. 

(Embedded business configuration of a computer) next, embedded - business - the 
concrete configuration of a computer is explained, drawing 2 - embedded [ this ] - business 
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■ it is the outline block diagram having shown only the configuration related to the 
identification information embedded processing to subject-copy image data among the 
hardware of a computer, it is shown in drawing 2 - as embedded [ this ] - business the 
computer has CPU1 each other connected by Bus B, the disk unit 2 for an input, ROM3, 
and the disk unit 4 for an output. 

[0039] The disk units 2 for an input are the hard disk drive unit which inputs the 
subject-copy image data 21 into CPUl according to the directions from CPUl, a floppy disk 
drive unit, optical -magnetic disc equipment, etc. 

[0040] ROM3 as an embedded function maintenance means and a computer-readable 
medium is a read-only memory holding the identification information embedded program 
31 performed by CPUl, identification information 32, and the embedded function 33. 
[0041] CPUl - embedded - business - it is the processor which performs control of the 
whole computer, and by executing the identification information embedded program 31 
read from ROM3, the orthogonal transformation processing section 11, the multiplier 
extract processing section 12, the multiplier permutation section 13, the multiplier 
embedding part 14, and the reverse orthogonal transformation section 15 are developed to 
the interior, and identification information embedded processing which shows an outline in 
drawing 3 is performed. In addition, the dotted line in drawing 2 shows the data flow in 
CPUl. 

[0042] The orthogonal transformation processing section 11 as an orthogonal 
transformation means performs orthogonal transformation processing mentioned above to 
the subject-copy image data 21 of the NxM pixel read from the disk unit 2 for an input, and 
computes the weighting factor to the combination of the NxM individual of a basis function, 
respectively. The weighting factor of this NxM individual forms the matrix of a NxM 
individual as well as a subject-copy image. The matrix of this weighting factor is hereafter 
called "multiplier distribution image data" for convenience. The orthogonal transformation 
processing section 11 notifies this multiplier distribution image data to the multiplier 
extract processing section 12 and the multiplier embedding part 14, respectively. 
[0043] The multiplier extract processing section 12 extracts the weighting factor of L pieces 
which should embed each signal of identification information from the multiplier 
distribution image data notified from the orthogonal transformation processing section 11, 
and notifies it to the multiplier permutation section 13. 

[0044] The multiplier permutation section 13 as a multiplier permutation means reads 
identification information 32 and the embedded function 33 from ROM3. And for every 
signal of identification information 32, the embedded function 33 is counted backward and 
all the input values of the embedded function 33 which makes the value of the signal an 
output value are calculated. And out of the calculated input value, a difference with the 
value of the corresponding weighting factor notified from the multiplier extract processing 
section 12 specifies what becomes the smallest, and notifies to the multiplier embedding 
part 13 as a rewriting value of the multiplier concerned. 

[0045] The multiplier embedding part 14 as a multiplier permutation means overwrites 
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the value of the weighting factor corresponding to the rewriting value notified from the 
multiplier permutation section 13 in the multiplier distribution image data received from 
the orthogonal transformation processing section 11 with this rewriting value. 
[0046] The reverse orthogonal transformation section 15 as a reverse **** conversion 
means performs reverse orthogonal transformation to the weighting factor of the NxM 
individual received from the multiplier embedding part 14. This reverse orthogonal 
transformation section 15 calculates total of a product with the value of the weighting 
factor corresponding to class doubling the basis function used in the orthogonal 
transformation processing section 11, and each basis function to the location of that pixel 
for every pixel of an identification information embedded settled image, and calculates the 
brightness value of that pixel. The reverse orthogonal transformation section 15 sends out 
the identification information embedded settled image 51 which consists of each pixel 
which did in this way and asked for brightness to the disk unit 4 for an output. 
[0047] The disk units 4 for an output are the hard disk drive unit with which the 
identification information embedded settled image data 41 passed from CPU1 is written in, 
a floppy disk drive unit, optical -magnetic disc equipment, etc. 

[0048] Drawing 3 is a flow chart which shows the contents of the identification information 
embedded processing performed by CPU1 which read the identification information 
embedded program 31. this identification information embedded processing - embedded - 
business - it starts ignited by an identification information embedded command being 
inputted through the keyboard which was connected to the computer and which is not 
illustrated. In S001 first performed after this identification information embedded 
processing starts, CPUl reads L signals Si [i=l"L] which make identification information 
32 from ROM3. 

[0049] In the following S002, CPUl reads L embedded functions 33 (fi(x) fi=l ■ L]) from 
ROM3. In the following S003, CPUl reads the subject-copy image data of a NxM pixel from 
the disk unit 2 for an input. 

[0050] In the following S004, CPUl performs orthogonal transformation processing 
mentioned above to the whole subject-copy image data read in S003, and generates the 
multiplier distribution image data of a NxM dot. 

[0051] In the following S005, CPUl chooses the weighting factor Ci [i=lL] in the location 
corresponding to each signal Si of the identification information read in S001 [i=l-L] from 
each weighting factors which constitute the multiplier distribution image data generated 
in S004. 

[0052] In the following S006, CPUl initializes the variable i for specifying the weighting 
factor Ci of a processing object, and sets it to "1." Next, CPUl performs loop -formation 
processing of S007 thru/or SO 10. Entering this loop formation, by S007 of the beginning, 
CPUl calculates the solution cik (l <=k<=Ki) of Ki individual which fills the relation of 
fi(cik) =Si based on Signal Si and the embedded function fi (x) of the identification 
information corresponding to the weighting factor Ci of a processing object. It is defined as 
the number of the solution cik of embedded function fi(x) =Si prepared here corresponding 
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to the i-th above-mentioned signal Si being Ki. 

[0053] In the following S008, although CPU1 is the closest to the weighting factor Ci 
concerned among each the S solution cik asked for the value in S007, it rewrites to a value 
the weighting factor Ci of the processing object in the multiplier distribution image data 
generated in 004. 

[0054] In the following S009, CPU1 confirms whether processing of S007 and S008 was 
performed to all the signals that constitute whether Variable i reached L and identification 
information. And if Variable i has not yet reached L, after incrementing Variable i in SO 10, 
processing is returned to S007. 

[0055] On the other hand, if Variable i has reached L, in SOU, CPU1 will perform reverse 
orthogonal transformation mentioned above to the multiplier distribution image of S008 
with which rewriting was made, and will generate the identification information embedded 
settled image data of a NxM dot. 

[0056] In the following SO 12, CPU1 writes the identification information embedded settled 
image data generated in SOU in the disk unit 4 for an output. By the above, CPU1 ends 
embedded processing. 

(Configuration of the computer for an extract) Next, the concrete configuration of the 
computer for an extract is explained. Drawing 4 is the outline block diagram having shown 
only the configuration related to the identification information extract processing from 
identification information embedded settled image data among the hardware of this 
computer for an extract. As shown in drawing 4 , this computer for an extract has CPU1 
each other connected by Bus B, the disk unit 2 for an input, ROM3, and the output unit 5. 
namely, the hardware configuration of the computer for an extract - embedded - business 
it is the same as that of the hardware configuration of a computer, the point that the 
identification information extract program 34 is stored in ** ROM 3, and the disk unit 4 for 
an output are not indispensable, and only the points that an output unit 5 is indispensable 
differ, therefore - if the disk unit 4 for an output and the output unit 5 are connected to 
Bus B while the both sides of the identification information embedded program 31 and the 
identification information extract program 34 are stored in ROMS - one computer - 
embedded business - it can be made to function as a computer and a computer for an 
extract 

[0057] In drawing 4 , the disk units 2 for an input are the hard disk drive unit which inputs 
the processing-object image data (identification information embedded settled image data) 
22 into CPU1 according to the directions from CPU1, a floppy disk drive unit, 
optical-magnetic disc equipment, etc. 

[0058] ROM3 as an embedded function maintenance means and a computer~readable 
medium is a read only memory holding the identification information extract program 34 
performed by CPUl, and the embedded function 32. this embedded function 32 - 
embedded business - it is completely the same as that of the thing of a computer. 
[0059] CPU1 is a processor which performs control of the whole computer for an extract, by 
executing the identification information extract program 34 read from ROM3, develops the 
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orthogonal transformation processing section 16, the multiplier extract processing section 
17, and the identification information calculation section 18 to the interior, and performs 
identification information extract processing which shows an outline in drawing 5 . In 
addition, the dotted line in drawing 4 shows the data flow in CPU1. 

[0060] The orthogonal transformation processing section 16 as an orthogonal 
transformation means performs orthogonal transformation processing mentioned above to 
the processing-object image data 22 of the NxM pixel read from the disk unit 2 for an input, 
and computes the weighting factor to the combination of the NxM individual of a basis 
function, respectively. And the multiplier distribution image data which consists of a 
weighting factor of a NxM individual is notified to the multiplier extract processing section 
17. 

[0061] The multiplier extract processing section 17 as a fetch means extracts the weighting 
factor of L pieces where each signal of identification information may be embedded from 
the multiplier distribution image data notified from the orthogonal transformation 
processing section 16, and notifies it to the identification information calculation section 18. 
[0062] The identification information calculation section 18 as a calculation means reads 
the embedded function 33 from ROM3. And the output value of the embedded function 33 
is calculated for each [ which was notified from the multiplier extract processing section 
17 ] weighting factor of every. And each calculated output value is arranged according to 
an array in the multiplier distribution image data of the weighting factor corresponding to 
the output value, and it outputs to an output unit 5. 

[0063] An output unit 5 is a printer which prints the display unit or output value which 
displays the output value of L pieces passed from CPUl. Drawing 5 is a flow chart which 
shows the contents of the identification information extract processing performed by CPUl 
which read the identification information extract program 34. This identification 
information extract processing is started ignited by an identification information global 
command being inputted through the keyboard which was connected to the computer for 
an extract and which is not illustrated. 

[0064] In S101 first performed after this identification information extract processing 

starts, CPUl reads L embedded functions 33 (fi(x) [i=l - L]) from ROM3. 

[0065] In the following S102, CPUl reads the processing-object image data of a NxM pixel 

from the disk unit 2 for an input. In the following S103, CPUl performs orthogonal 

transformation processing mentioned above to the whole processing-object image data read 

in S102, and generates the multiplier distribution image data of a NxM dot. 

[0066] In the following S104, CPUl chooses the weighting factor Ci [i=l"L] in the location 

corresponding to each signal Si of the identification information 32 in an embedded 

computer [i=l-L] from each weighting factors which constitute the multiplier distribution 

image data generated in S103. 

[0067] In the following S105, CPUl initializes the variable i for specifying the weighting 
factor Ci of a processing object, and sets it to "1." Next, CPUl performs loop -formation 
processing of S106 thru/or S108, Entering this loop formation, by S106 of the beginning, 
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CPU1 calculates the output value Si over Ci of the embedded function fi corresponding to 
the weighting factor Ci of a processing object (x). 

[0068] In the following S107, CPU1 confirms whether processing of S 106 was performed to 
all the weighting factors where whether Variable's i having reached L and identification 
information may be embedded. And if Variable i has not yet reached L, after incrementing 
Variable i in S108, processing is returned to S106. 

[0069] On the other hand, if Variable i has reached L, in S109, CPU1 will arrange in as an 
array within the multiplier distribution image data of each corresponding weighting factor 
Ci [i=l-L] all the output values Si [i=l-L] calculated in S106, and will output them to an 
output unit 5. Consequently, an output unit 5 can display or print the data corresponding 
to identification information, if the processing-object image data read in S102 is 
identification information embedded settled image data. 

[0070] In addition, as orthogonal transformation used for this operation gestalt, a 
two-dimensional discrete cosine transform (two-dimensional DCT), two-dimensional 
dispersion sign conversion (two-dimensional DST), or a two-dimensional Hadamard 
transform is employable. Next, the contents of concrete processing of the identification 
information embedded processing at the time of adopting two-dimensional DCT as 
orthogonal transformation and identification information extract processing are explained 
below as an example 1. 
[0071] 

[Example l] In this example 1, subject-copy image data presuppose that it consists of NxN 
(however, N»8) pixels, as shown in drawing 10 (a). Moreover, the brightness value of each 
pixel of subject-copy image data presupposes that they are the gray scale of 0-255. 
Moreover, identification information presupposes that it is the image data (henceforth 
"signature image data") which gave the white brightness value (255) or the black 
brightness value (0) alternatively to each 8x8 pixel, and expressed the 
European-languages character "FJ", as shown in drawing 10 (d). In connection with this, 
identification information embedded settled image data shall be called "signature image 
embedded settled image data." Moreover, since the weighting factor obtained by 
two-dimensional DCT is equivalent to the reinforcement (amplitude) of each frequency 
component in subject-copy image data, it shall call for convenience the "multiplier 
distribution image data" mentioned above "frequency distribution image data." Moreover, 
the periodic function which followed serrate [ which is shown, function f (x), i.e., the graph 
of drawing 10 (e), of one ** applied in common to each weighting factor as an embedded 
function, ] shall be prepared. 

(Identification information embedded processing) Drawing 6 and drawing 8 are flow charts 
which show the identification information embedded processing by the example 1. 
[0072] In S201 of drawing 6 first performed after this identification information embedded 
processing starts, CPU1 reads signature image data S (i, j) [i=0-7, j=0-7] from ROM3. 
[0073] In the following S202, CPU1 reads embedded function f (x) from ROM3. In the 
following S203, CPU1 reads the subject-copy image data of a NxN pixel from the disk unit 
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[0074] In the following S204, CPU1 gives two-dimensional DCT to the whole subject-copy 
image data read in S203, and generates the frequency distribution image data of a NxN 
dot as shown in drawing 10 (b). Specifically by these S204, CPU1 performs the 
two-dimensional DCT processing subroutine shown in drawing 7 . 

[0075] It goes into this two-dimensional DCT processing subroutine, and in S301 of the 
beginning, CPUl initializes the variable i which shows the train (let 0 be the leftmost) in 
the inside of the frequency distribution image data of weighting-factor C for calculation (i, 
j), and sets it to "0." 

[0076] In the following S302, CPUl is set up as the function Wi of Variable x (x) is shown 
in the following type (2). Variable x supports the train (let 0 be the leftmost) of each pixel in 
subject-copy image data here. 
[0077] 

[Equation 2] 

[0078] In the following S303, CPUl initializes the variable j which shows the line (0 is 
made into the best) in the inside of the frequency distribution image data of 
weighting-factor C for calculation (i, j), and sets it to "0." In the following S304, CPUl is set 
up as the function Wj of Variable y (y) is shown in the following type (3). Variable y 
supports the line (0 is made into the best) of each pixel in subject-copy image data here. 
[0079] 

[Equation 3] 

[0080] In the following S305, CPUl performs the following type (4) based on the brightness 
value G (x y) of each pixel contained in both the functions Wi (x) set up in this time, and Wj 
(y) and subject-copy image data, and makes the computed series the value of 
weighting-factor C (i, j) specified by the current variables i and j. 
[0081] 

[Equation 4] 

[0082] In the following S306 ? CPUl confirms whether the value of the variable j at present 
amounts to (N-l). And when the value of Variable j does not yet amount to (N-l), CPUl 
increments Variable j in S307. Then, in S308, CPUl substitutes the variable j in this time 
for the following type (5), and resets it up as a new function Wj of Variable y (y). 
[0083] 

[Equation 5] 
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[0084] CPUl is after that, returns processing to S305 and computes weighting-factor C (i, 
j) which consists in the following line. On the other hand, when it judges with the value of 
the variable j at present amounting to (N-l) in S306, CPUl advances processing to S309. 
In these S309, CPUl confirms whether the value of the variable i at present amounts to 
(N-l). And when the value of Variable i does not yet amount to (N-l), CPUl increments 
Variable i in S310. Then, in S311, CPUl substitutes the variable i in this time for the 
following type (6), and resets it up as a new function Wi of Variable x (x). 
[0085] 

[Equation 6] 

[0086] CPUl is after that, returns processing to S303 and computes weighting-factor C (i, 
j) which consists in the following train. On the other hand, when it judges with the value of 
the variable i at present amounting to (N-l) in S309, CPUl judges that all the weighting 
factors that constitute frequency distribution image data were computed, ends this 
subroutine, and returns processing to the main routine of drawing 6 . 

[0087] In the main routine of drawing 6 by which processing was returned, processing is 
advanced to S205 after completion of S204. The processing to S212 is after these S205 
processing for taking out the field (field of i=0"7 and j=0~7) ( drawing 10 (c)) which shows 
the field and bottom frequency component which are shown in drawing 10 (b) thru/or (g) ? 
and which show DC component in frequency distribution image data ( drawing 10 (b)) like, 
and embedding signature image data to this field. Thus, the effect to which having 
embedded signature image data only at the weighting factor which shows the weighting 
factor which shows DC component, and a low-frequency component gives fluctuation of 
some of DC components and low frequency components to the image quality of signature 
image embedded settled image data is because it is small. 

[0088] In S205, CPUl initializes the variable i which shows the train in the inside of the 
frequency distribution image data of weighting-factor C for fetch (i, j), and the train of the 
pixel S for reference in a signature image (i, j), and sets it to "0." 

[0089] In the following S206, CPUl initializes the variable j which shows the line in the 
inside of the frequency distribution image data of weighting-factor C for fetch (i, j), and the 
line of the pixel S for reference in a signature image (i, j), and sets it to "0." 
[0090] the following S207 - CPUl - the signature image data ( drawing 10 (d)) from ROM3, 
and embedded function f (x) - (- drawing 10 (e)) is read and the total input value cijk of 
embedded function f (x) which makes an output value the brightness value of the pixel S 
for reference in the signature image data specified by the variable i at present and 
Variable j (i, j) is calculated. That is, the solution cijk (l <=k<=Ki) of the Kij individual 
which fills fvcijk) =S (i, j) is calculated altogether. Here, the number of Solution cijk to 
f(cijk) =S (i, j) is defined as a Kij individual. 

[0091] In the following S208, CPUl chooses the thing nearest to the value of 
weighting-factor C for fetch (i, j) specified by the variable i at present and Variable j among 
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the total input values cijk calculated in S207. And the value of weighting-factor C for fetch 
(i, j) specified by the variable i at present and Variable j is transposed to the selected input 
value cijk (refer to drawing 10 (£)). 

[0092] In the following S209, CPU1 confirms whether the value of the variable j at present 
amounts to "7." And when the value of Variable j does not yet amount to "7", CPU1 returns 
processing to S207, in order to permute weighting-factor C (i, j) of the following line, after 
incrementing Variable j in S210. 

[0093] On the other hand, when it judges with the value of the variable j at present 
amounting to "7" in S209, CPU1 advances processing to S211. In these S211, CPU1 
confirms whether the value of the variable i at present amounts to "7." And when the value 
of Variable i does not yet amount to "7", CPU1 returns processing to S206, in order to 
permute weighting-factor C (i, j) of the following train, after incrementing Variable i in 
S212. 

[0094] On the other hand, when it judges with the value of the variable i at present 
amounting to "7" in S211, CPU1 advances processing to S213. In these S2 13, CPUl gives 
two-dimensional reverse DCT to the whole frequency distribution image data ( drawing 10 
(g)) containing weighting-factor C (i, j) by which the value was permuted in S208, and 
generates the signature image embedded settled image data of a NxN dot as shown in 
drawing 10 (h). Specifically by these S213, CPU1 performs the two-dimensional reverse 
DCT processing subroutine shown in drawing 8 . 

[0095] Going into this two-dimensional reverse DCT processing subroutine, CPU1 defines 
the function used for this processing by S401 of the beginning as following type (7) - (10). 
[0096] 

[Equation 7] 

[0097] In the following S402, CPUl initializes the variable x which shows the train (let 0 
be the leftmost) within the signature image embedded settled image data of the pixel R for 
calculation (x y), and sets it to "0." 

[0098] In the following S403, CPUl initializes the variable y which shows the line (0 is 
made into the best) within the signature image embedded settled image data of the pixel R 
for calculation (x y), and sets it to "0." 

[0099] In the following S404, CPUl computes the brightness value of the pixel R of the 
signature image embedded settled image data specified by the variable x in this time, and 
y (x y). Specifically, CPUl is substituted for the function of the formula (10) which defined 
the value of the variable y in this time by S401 while it assigns the value of the variable x 
in this time to the function of the formula (8) which S401 defined. Based on the function of 
the formula (7) defined by S401 as the function of a formula (8) with which CPUl moreover 
performed the above-mentioned substitution and the function of a formula (10), and the 
list, and the function of a formula (9), the following type (ll) is performed and let the 
computed series be the brightness value of the pixel R (x y) specified by the current 
variable x and y. 
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[0100] 

[Equation 8] 

[0101] In the following S405, CPU1 confirms whether the value of the variable y at present 
amounts to (N-l). And when the value of Variable y does not yet amount to (N-l), after 
CPU1 increments Variable y in S406, it returns processing to S404 and computes the 
brightness value of the pixel R (x y) which consists in the following line. 
[0102] On the other hand, when it judges with the value of the variable y at present 
amounting to (N-l) in S405, CPU1 advances processing to S407. In these S407, CPU1 
confirms whether the value of the variable x at present amounts to (N-l). And when the 
value of Variable x does not yet amount to (N-l), after CPUl increments Variable x in S408, 
it returns processing to S403 and computes the brightness value of the pixel R (x y) which 
consists in the following train. 

[0103] On the other hand, when it judges with the value of the variable x at present 
amounting to (N-l) in S407, CPUl judges that the brightness value of all the pixels that 
constitute signature image embedded settled image data was computed, ends this 
subroutine, and returns processing to the main routine of drawing 6 . 

[0104] In the main routine of drawing 6 by which processing was returned, processing is 
advanced to S214 after completion of S213. In these S214, CPUl outputs the signature 
image embedded settled image data of a NxN pixel to the disk unit 4 for an output. 
(Identification information extract processing) Drawing 9 is a flow chart which shows the 
identification information extract processing by the example 1. 

[0105] In S501 of drawing 9 first performed after this identification information embedded 
processing starts, CPUl reads embedded function f (x) from ROM3. In the following S502, 
CPUl reads the processing- object image data of a NxN pixel as shown in drawing 11 (a) 
from the disk unit 2 for an input. 

[0106] In the following S503, CPUl gives two-dimensional DCT to the whole 
processing-object image data read in S502, and generates the frequency distribution image 
data of a NxN dot as shown in drawing 11 (b). Specifically by these S503, CPUl performs 
the two-dimensional DCT processing subroutine shown in drawing 7 . 

[0107] The processing to S510 is after the following S504 processing for taking out the field 
(field of i=0-7 and j=0*7) ( drawing 11 (c)) which shows the field which is shown in drawing 
11 (b) thru/or (e), and which shows DC component in frequency distribution image data 
( drawing 11 (b)) like, and a low frequency component, and extracting signature image data 
from this field. 

[0108] In S504, CPUl initializes the variable i which shows the train in the inside of the 
frequency distribution image data of weighting-factor C for fetch (i, j), and sets it to "0." In 
the following S505, CPUl initializes the variable j which shows the line in the inside of the 
frequency distribution image data of weighting-factor C for fetch (i, j), and sets it to "0." 
[0109] In the following S506, it computes eye i train and output-value S' (i, j) of the j-th line 
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by CPU1 reading an embedded function (the same as that of drawing 11 (d) and drawing 
10 (e)) from ROM3, and assigning the value of weighting-factor C for an extract (i, j) 
specified by the variable i at present and Variable j to embedded function f (x). 
[0110] In the following S507, CPU1 confirms whether the value of the variable j at present 
amounts to "7." And when the value of Variable j does not yet amount to "7", CPU1 returns 
processing to S506, in order to compute output-value S' of the following line, after 
incrementing Variable j in S508. 

[01 11] On the other hand, when it judges with the value of the variable j at present 
amounting to "7" in S507, CPU1 advances processing to S509. In these S509, CPUl 
confirms whether the value of the variable i at present amounts to "7." And when the value 
of Variable i does not yet amount to "7", CPUl returns processing to S505, in order to 
compute output- value S' (i, j) of the following train, after incrementing Variable i in S510. 
[0112] On the other hand, when it judges with the value of the variable i at present 
amounting to "7" in S509, CPUl advances processing to S511. In these S511, CPUl 
outputs the image data which consists of an output value (brightness value) of the ixj dot 
computed in S506 to an output unit 5. If processing-object image data is signature image 
embedded settled image data at this time, this image data will turn into the same image 
data ( drawing 11 (e)) as a signature image ( drawing 10 (d)). 
[0113] 

[Effect of the Invention] According to this invention constituted as mentioned above, 
without worsening the image quality of image data not much, from a third person, it is the 
gestalt which is not discriminable and identification information can be embedded to 
subject-copy image data. And the embedded identification information can be extracted 
without subject-copy image data. Therefore, since the rightful claimant of image data etc. 
does not need to save and manage subject-copy image data, he does not need to prepare 
mass storage. 



DESCRIPTION OF DRAWINGS 



[Brief Description of the Drawings] 

[Drawing ll The principle Fig. of this invention 

[Drawing 2] embedded [ by the 1st operation gestalt of this invention ] - business - the 
block diagram showing the outline configuration of a computer 

[Drawing 3] The flow chart which shows the contents of the identification information 
embedded processing performed by CPUl of drawing 2 

[Drawing 4] The block diagram showing the outline configuration of the computer for an 
extract by the 1st operation gestalt of this invention 

[Drawing 5] The flow chart which shows the contents of the identification information 
extract processing performed by CPUl of drawing 4 

[Drawing 6] The flow chart which shows the contents of the identification information 
embedded processing by the example 1 
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[Drawing 7] The flow chart which shows the contents of the two-dimensional DCT 
processing subroutine performed in S204 of drawing 6 

[Drawing 8] The flow chart which shows the contents of the two-dimensional reverse DCT 
processing subroutine performed in S213 of drawing 6 

[Drawing 9] The flow chart which shows the contents of the identification information 
extract processing by the example 1 

[Drawing 101 The explanatory view showing the flow of the signature image embedded 
processing by the example 1 

[Drawing 111 The explanatory view showing the flow of the signature image extract 
processing by the example 1 
[Description of Notations] 

1 CPU 

2 Disk Unit for Input 

3 ROM 

4 Disk Unit for Output 

5 Output Unit 

11 Orthogonal Transformation Processing Section 

12 Multiplier Extract Processing Section 

13 Multiplier Permutation Section 

14 Multiplier Embedding Part 

15 Reverse Orthogonal Transformation Processing Section 

16 Orthogonal Transformation Processing Section 

17 Multiplier Extract Processing Section 

18 Identification Information Calculation Section 
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O^TOAtffS^**^o fit, >K*?>n/cA*fiO 

^e>, ^mfflm^js^i 2^e>aaisnrc*fis-rss 

[0045] ««««i^ai: LTo«»aii* 1 4 
*o«»a«8Pi 3*^5afti*nfe»»tfjc:*rjE'rsa 
[0046] m%w£W£m£hT<Dmm£&wmi 5 

5fcJ\ MgU1f«aiii5SiS«©S-H*@^, 
S^gp 1 5 te\ c CD J; *5 tc bTJi^^fefc^B*^ 

6a««gijiWBffiiiWH«5 1*, ai^fvx^gi 

4tc^Dm~To 

[0 0 4 7] tMfflfVX^gl4^ CPU 1^58 

zftrcmmmmmj&mwMy'-z 4 i^#^jA^na^ 

[0 0 4 8] H3B\ M9J1»«aii^n^A3 1 «»! 

ft&k,itc p u i ^oT^sns^gimwffiiiffiia 



15 



(9) 



f# P*T¥ 10-234012 

16 



h Lrc&§kfflfc£ff*iEti% s o o i ra, cpunt 

MS'J1f*3 2*&"TLffl(Dm#Si [i = l~L] £\ 
ROM3^ei^o 

[0 0 4 9] &0)S 0 0 2T«, CPU Hi, LfflO±l 
JAMtS(3 3 (f t(x) [i = 1-L] ) *ROM3i^R 
&7&tS 0 %L<D S 0 0 3 CPUia, NX MII(D 

iSHflfrr— A^ffl^Vx^g2;^5B^3&t? 0 

[0 0 5 0] $1<DS 0 0 4T&, CPUlli, S00 3 

mmm&m it, n x m f v y <^%mmm^rf— $ 

[0 0 5 1] ^OSOOSm CPU Hi. S004 

s o o i \cxmwhh;rm®m®.<o&m^ 
Si [ i = i ~l: tc*f/£-rsffiatcssa*«»c 

x [i - l-L] *8S?-r5 0 

[0052] ^soo6 cpunt ffla*f^ 

T "1" ^t§o ^tc s CPU Hi, S 00 7MS0 
tD S 0 0 7 C P U Hi, 98mtm<OWfrffi%L C i 

*3vt, f i ( c ik ) = s i (DmmztmrcT k ± m^m c * 

(l^k^KO ^S<a6§o ClIC ±fdifg©i§ 
S i fc^jSLTffljftStifcSBiHa f i (x) - S i(DM c ik 

[0 0 5 3] ^S 0 0 8m CPUlti S00 4 

«»Ci*ffi*\ S 0 0 7fcTjR»£ftfc&j&cik *g 

[0 0 5 4] S 0 0 9 -VIZ, CPUUt i ^ 

^c^lts 007MS00 8<ojaig*HfTLrc5^g 

^*fx'^T8„ fit, *Sfti^LrlLT^ 

ftttntf. soiot:tig[i^>^u^>n/ci 
t\ jas* s o 0 7 tcM1~o 

[0 0 5 5] cmcMLXmk i ^Uu}gLT^tUi\ 
S011 fcSs^T. S 0 0 8 TCD»^cO* 

[0 0 5 6] $i<D S 0 1 2 C P U Hi. S 0 1 1 



3>tfa.-^t±, mnc^xBie^oTg^^nrcc p 

Ul, AifflrVX^gf2, ROM 3, Rtfffi^jSK 

ra-i?* d , *r o m 3 fyicmmmmm&^u 3 

10 ^T^<ffi±f^g5^^T^^^CD^^S^§ 0 fi£ 
oT. ROM3rttCi|glJfl|«Sii^ , n^A3 lfitfl 

[0 0 5 7] S4ttVT, AAfflfVX^ii2^ 
C P U 1 ^StDJB^fcjSCTfflaj^Bli^-^ (WiiJ 

mmm&ffim^-z) 22%cpui ica^-ts^— 

X^B, fpftSo 

[0058] ^H»«^sat;3>if a-^prii« 

#^LTOROM3^ C P U UcTHfrSftSiBgfJIi 
?gaffi7 P n^A3 4, StfSEil|»3 2*«J*LTV> 

5S*(BLSffl^^ut?$5o c<d«3&b§»3 2&. m 

[0 0 5 9] CPU Hi. mmm^W*.-#^#<Dffl 
rclffiS'Jt«Wm^n^A3 4*Hfr"rscfcfc:<fcr? 

30 r. *<Dftmcm£g&%mgti 6, «»fflffifflaa5i 
■mBU««ttajaa**ffrso a*. B4fc:*5tf3j& 

*8«x C P U 1 fi^fe^§f-3r«n^LTl/^o 
[0 0 6 0] B0?g*^Sfc LT<Dfi3^g|ffia^ 1 6 
fr±* A*ffl-r>TX^^B2*^K*HlLfeNXMH*<D 
liSftlif-^ 2 2 tc3* LT±JELfcit^lSlffla 
^ffLT, Sj£l8iSC£DNxMfflOffi^tC^-r^a^ 

40 T-^o 

[0 0 6 1] ^m¥S^ LTO«»»tB«ra» 1 7 tt. 

it^^a^ 1 6 t>m9®-ztircffi?&ftif;m»T—* 

5LtOi»»»«»jULT, KBiJfiBBJMJSP 1 8-\il 

[0 0 6 2] Stb#©^LTOSSWIt«*ffig!3l 8 
R0M3^6aaBBiR3 3*R&aty 0 ^LT, «iStft 
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[0 0 6 3] ffiMISii, CPUlfr&SStlTcLffl 
A3 4^rM^?AA,rcC PU 1 tc£ ^T^f^tlSMJSUli 

wcm-j s i o i -e&, cpuiit LfflcoajA 

88^3 3 (f i(x) [i = 1— L] ) *ROM3fr5i* 
[0 0 6 5]*OS102m CPUliis NXMI 
tc^iifro ^£0 S 1 0 3 CPUlii, S 1 0 2 

-**4filtt5o 20 
[0 0 6 6] &<7)S 1 0 CPU Hi, SI 03 

^(S^Si [ i = 1 ~L] tc^/£-rs{ftStcfe§S^« 

8!(Ci [i = i~l] *as?fSo 

[0 0 6 7]^(DS105 CPUHi, jffiffiK* 

T " 1" i:-r^ 0 CPU Hi, S106MS1 

<d s i o 6 c p u is. mmmcDWLfrmiiCi 30 

[006 8]»S10 7^fci\ C P U Hi, ^ i ^ 

^ttOfeS^T^M^StC^fbTS 1 0 6^3* 

Uc^gLTV^ttUf, S 1 0 81CT^8S[ i ^M'^U 
> h Lfcf£T\ MJl^r S 1 0 6 tcMTo 

[0 0 6 9] cnicfthrmmitfLicmLT^tii^ 

CPU Hi, S 1 0 9 tcfel^T, S106tcT*46 5tL 40 

fc£T<Dm?jffiSi [ i = i~l] £\ ^ter^^a^ 
MCi [ i = i ~l] omm^mmy'-zftvcDm 

s»5«. s i o zicrm^^rmmMmmm'r-^ 

[0 0 7 0] **«»JB^fflV^n5iB3!^lRfc 
LTtt. Z^TcJBftn^vS* (2^DCT) , 2 



[0 0 7 1] 

[§»#(! 1] cdtiiim HMft-r— 01 
0 (a) ic^T&oK, NXN (fSU N>>8) Bffi 

ffffifBS. 0-2 5 5£D^V-X^r— ;l/t$5tt8 0 
Sfe, WSU«*fci\ 110 (d) 8X8 

fflcD«-B*tcaoiififii (255) xai©wai 

(0) *S&ift£##LTHfcSC*: Tf Jj ^SLfcBff 

DCTt±oTS6h5i»«»t mm»T~-#\c& 

(x) , Hi 0 (e) (D^^yic^ti^mm^R 
[0072] c (Dmmmmm&mmtfx ^-m mm 

W£*fTSttSH6<OS 2 0 1 -Ctt. CPUlti i4 
H«T^£S(i,j) [ i -0 — 7, j - 0 — 7] £\ R 

[0 0 7 3]*(OS 2 0 2T1i, C P U H3\ ffiAMifc 
f (x)*R0M3fr&K#i&tS' o #OTS 2 0 3T^, C 
PUIS, NxNIiOglir-^^ AJjJBtVX 
^g2^£MW&t? 0 

[0 0 7 4]^S 2 0 4m CPU Hi, S203 

icTm&&/vtmmMT—$±mctthT 2ykjiD c t 

^rSSLT. 810 (b) ^*ti4NxN Vy h(DMtifr 
4T?«. CPU Hi, 0 7tC^-T2^DC Tffiil^^ 
[0 0 7 5] C0 2^DCT®W7;^f>lcA-3 

rgfflo s 3 o i c p u Hi, *ffi*f*s*«a 

c a jjosir^itf-^ ^t^ij (o^a^t 
•rs) *^-r^Ri*. sifflffcLr "o" t-rso 

[0 0 7 6] ^£DS 3 0 2 Tti, C P U Hi, ||i(xO 

itWi(x)^ ties (2) ^-raosss-rso c 

[0 0 7 7] 
KSC2] 

Wc(x)= -/J (2) 

[0 0 7 8] &<DS 3 0 3 m C P U Hi, *m*fft 
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*«±t*r«) *sts»j*. smmllt "o" £-r 

§o MS 3 0 4 CPUl^ My OIB^W 

j(y)*> TIES (3) t^TIDiSt^o CCtf 

[0 0 7 9] 

1 

C(i,j)= I z G(x,y)-Wi(x 

x=0 y-Q 



[0 0 8 0]»S 30 5m CPU1B, ?P»£;tC 
fc^T8^jntl^BH»Wi(x), Wj(y)^J3giffi« 
X"-^tc#sns#H*0»S#G (x,y) fctcg^T 

TtSxt (4) *nffL. KaisnfeiisaR*. uteris 

[0 0 8 1] 
[«4] 



•Wj(y) (4) 



[0 0 8 2] 3 0 6^ CPU1& 

xym e fit, *7£^»j Ott^(N- DtcML 
TV^^i^C^ C P U Hi. S 3 0 7 tC&'l^T. l£ 
iSCj *^>^U^>h"r«o Sc^T. CPU Hi, S3 
0 8t*t^ *TSBa£ (5) fcftA 

U ffiRyO*ffe4H»W J (y)4:LTaa^Lie'ro 20 

[0 0 8 3] 

[«5] 

[0 0 8 4] C PU 1 H\ *£M£^ ffiil^S 3 0 5fC 
ML. ^cofftc#TSS*«fRC(i,j)*MEti}1-«o c 

Wi 00 = cos i 

[0 0 8 6] CPUUt fflS^ S 3 0 3 ic 

ml. ^^jtc#-r5s*«ac(ij)**mrso £ 30 

LT^SfcS 3 0 9JCT3pJSLfc«^fc:tt. CPU1 

[0087] mmtfm^ftrcm6<D^s{>}i<^yiz$5 

V^T&. S 2 0 4£D^T^. fflStfS 2 0 5^ili6^n 
£ 0 lOS 2 0 2 1 2 0 1 0 

(b) nm (g) ^^^n^s^ N jss&^h*^- 

? (H10 (b)) ^ODCjS^^-TfliJKatficJSS 40 
»/£#*^-r®« ( i - 0 - 7_iio j = 0 ~ 7 ^M^) 
(HI 0 (c) ) ^DffiLt, caSWcSSBHK-r 

[0 088] S205 Vte, C P U Hi, 3Xffi?f^fi^ 

«»c ft j)ois»»ffiBiT-i"fi'eow, atf. 



LTV^Si: S 3 0 6tcT¥iMLfe«^c«. CPU1 
fcJU ffl-S* S 3 0 9 tC}li6&o CCD S 3 0 9 C P 

U 1 SiBfrStCfcttS^Sfc i <Offitf(N- 1 )tc3SbT 
^5*Jf5Mfxy^n o ^LT, *f£S» i £>ffi 
tf(N- l)tc^LTV^V^*^iia, CPU Hi, S3 

cpuiti, s 3 i i tcfev^r, ^jSTto^a i *r 

IBS C6) icftAU SixOi/tWWi(x)^L 
TS^ Lit "To 

[0 0 8 5] 

[»6] 

..... (6) 

SJSMfcLT "0" tT%o 
[0 0 8 9] &<DS 2 0 6m C P U Hi. ffiffifctft 

c a j ) <omwi»»^m»r t - $ *-e<z>fi\ r 
j mmitLT "o" t« 0 

[0 0 9 0] aoS 2 0 7T*B, C P U Hi, ROM3 
^&*&B«T f — # (Ell 0(d)) StfffiAHIfcf 
(x) (Ell 0(e)) Srft^ffiLT. a^tcfc^S^ 

#Sg*MfcHr*S (ij) C»®MB^m*ffli:-r^ffljA^ 
Sf (x) O^A^fScxjk ^*i6^ 0 f 

(cijk)=S (i, j) %mrc?Kij ffltO^Cijk (1 
^k^Ki) £\ ^TM&b%o clclT\ f (cijk ) = S 

(i, j ) icttTzmcij* (ommzKv mt^m^n 

[0 0 9 1 ] %i<D S 2 0 8 T-«, CPU1S, S 2 0 7 

c(ij)offifc:st,affv^fc©«ai(? , r5o *lt, sp# 

jSfctstf«£ft i RtfSSc j »cj:-3TWJ£2ft*8mj*f 
*fi*Hm C (i j) Offl*, BB?£nfcA*« c ijk EI 
#Sl?L£ (Ell 0(f )#H» o 



( 12 ) f# Pwl¥ 1 0-2 3 4 0 1 2 

21 22 

[0092]*OS 20 9T?tt, C P U 18, SH?j6fc tt, C P U 1 tt, ^^TfOffiR x Oil* S401ICT 

*»S^BCJ OfflCW "7" taLTt/>S^if5fr*^x SBLfci? (8) ©MRteftATS 1 1 1 g!Btj£T* 

>y^-T§o f tt> *f£SSSCj Ot^ "7" t£iiLTV> ©^gfcyOfit&S 4 0 1 tcT^BLfciS (1 0) ©|I8£Sc 

fcVfc^Ktt, CPUltt, S 2 1 0fc*SV>Tg» j * fcftVTSo Z<D_tT\ CPUltt, ±fEftA£fTofc 

■YVfJ^yHftft-?, ^OfxOa*««C(ij)* 5£ (8) OigjRatfS (1 0) (ORft. Mt^C, S40 

B^-rsz-cfeic, san^s 2 o 7^\M-r= ncrmmhfc^ (7) ©wss^s o) omswcg 

[0 0 9 3] cm^LT, TIBxt (11) &£frU *W*nfciR»*, 

^ "7" tILT^5iS 2 0 9tTfiJ^Lfc^t ^&oaa»xatf y fcioTW£Sn«BRR (x,y) 

tt, c p u l tt, ©H* s 2 l l icmtbZo co) s 2 i ©ffffiffifc-f 5„ 

1 Ttt, CPUltt SP^JCfetfSSifc i io [0 10 0] 
"7" KlLTI/^^i^^x'^-r^o ^LT, B&8] 
**SSSRi ©ffitf "7" E2LTV^V^Ea, CP R(x,y)= I Z C (iJ)-Wi (x)-Wj(y) (11) 

U 1 tt, S 2 1 2 icfcl^T^lfc i^^'J^^hbft '"° '"° 

arc, *o?ijca*«ftc (i ,j)£g#rr§fci6fc, & [oioi]*os4os cruitt mmmc 

2 0 6-\M-To fcW-SffiRyOffltfCN- l)ClLn^^i:9)!)^f 

[0 0 9 4] cnt»Lt, SBfr&fcfcttSgft i ©ffi i^m. ?LT-, *Saaj[yOffiS«(N-l)li:2L 

"7" fc»LTV3 S 2 1 1 tTfiJ^L/£«^tC TVfcV^fctt, CPUltt, S 4 0 6 K*^TfR 

tt, C P U 1 tt, MS^rS 2 1 3s\m*b&o COS21 y^-T V^U^VhLfc&fc, %m%S 4 0 4^IU 

3 Ttt, CPUltt, S 2 0 8 tt«a<«j|?nftl* ^OfrKfiFf *B*R (x.y) v>nmmm-&t 
%&C(l,j)%^tsm®£lft1iiW»7 : —2 (Ml 0 20 [0 1 0 2] cntC*fLT, ^Bf^r&tfS^gryCDfii 

(g)) ^*fc*fLT2*7c5S*DCT*jSbT, El 1 0 ff(N-l)fcaLT^«fcS 4 0 5ETW^Lfci^t 

(h) tc^T«&NXN K-y h©a£B<£S3&^jg#f f tt, CPUltt, ffiU% S 4 0 7 ^\jt46§ 0 CCDS 40 

-££r£fiSc"rS 0 Mfr&^tt, COS 2 1 3 m CP 7Ttt, C P U 1 tt, ^P^Cfctf^iS! x ©fit*HN 

Ultt, ^8fc^t-2^^DCTffl@I+)-:7;l/— ^>«r - 1 ) IcjgLTl^Sfr t* 5 fr^i >y ^-TSo fLT, 

^WSo *7£^8txOffi#(N- 1 )fcaLT^ftV^§-g-tCtt, C 

[009 5] C©2^a5DCT^l-9-y;l/-f=->{cA P U 1 tt, S 408K8^T»[xS'f>'j"J^yH 

^>t*#j«s 4 o i -ett, cpuitt, cojast;:^ fc^t, mm&s a 03^11, wom^-r^mmR 
ztizmm*. ties (7) ~ (1 0) (D&vftmr ( x , y ) <DMmm*»tii?%o 

5 ° [0 10 3] CtUCfcfLT, ^HfjStCfeW-SggfcxOffi 

[0 0 9 6] 30 #(N-l)fc:jlLTV>5£S 4 0 7 {cTW^L/cif-g-fc 

CSt7 ^ a, c p u i tt, w&mmni&mwmT-t'tzMmTz, 

WoW=Vtt (7) £T©B*©«fitttf*msftfc£$iwfLT, 

W, 00- Vfr™ g N ' (8) 

k [0104] mmtfm-$nrcm6<D*^»v— ^>icts 

WoW-Vtt < 9 > i^Ttt, S 2 1 «yi*^S 2 1 4'\}ti6C>ft 

/I * (2y+1j «° COS2 14Ttt, CPUltt, NxNMJiOWS 

[ 0 0 9 7 ] &<o S 4 0 2 Ttt, CPUltt, 40 (aKSfflH^ffiftLSD H 9 tt, gffiffll 1 fc: ± 

S*R (x,y) OW^HiUI&MHfi^-^ftTO^J ftmjQa%St7n'-ft- hT?feS„ 

(OfcSEfTS) *jRTSHSx*, SJJtWfcLT "0" [0 10 5] ZoyWWffl&WMm&T.Z— h 

^" rs ° Mc*f?£tl£EI90S 5 0 1 Ttt, CPUltt, S3& 

[0098]^OS403 T?tt, CPUltt, Sffi*t#i Rft f (x)*R OM 3*»5«*iity 0 ^OS 5 0 2T 

HfSR (x >y ) com%mwitm&ffimm-T--2ftT°(Dff a, cpuia, @i 1 ( a ) icmT^v^Nxumm 

(ofcSLtfc-rs) *jRfg*y*, SBHfkLT "0" ©ft3*f*B«'r-**, Aifflf^X^gl2A^i 

[0 0 9 93 SOS 4 0 4 Trtt, CPUltt, SHS^T? [0106]*OS 503 Ttt, CPUltt, S 5 0 2 
f-^OI*R Ot.y) OBfii*IBt«. ^#WtC so DCT^ILT, 011 (b) tgfiSNxNF'yh 
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S 5 0 3 Tti, CPUl^ E( 7 tC^f 2 ^tcD C Tffi 

[0 10 7] &<DS 5 0 4KBS 5 1 0S^GDflEtfi«:. 
Ell 1 (b) TiS (e) fc^SftSfiifc:, 

(01 1 (b) ) ^DCl^tlfStf 
ffiH»*«»*^*T«« (i =0-7fiOj = 0-70 

mm) mi 1 c C ) ) aoffibt, co««^sa 

[0 10 8] S 5 0 4 C P U Hi. JfcaUti«^ 10 

«C (i *^<D5»J**-rS 
SJJHffcLT "0" £T%o *«S 5 0 5m 
C P U Hi, 3{tfl?Mfea*«»C (i j)*DHjft»»ffiH 
•■r-**^fT*^-rS»j*, SSJMbLT "0" 

[0109]^©S 5 0 6m CPUl^ ROM3 

frzw&MSi (mi i(d), eiio (e) tm—) &m 

(x) tcftALT. i^J@ • j ffSOtiitsms 9 (ij) 20 

[ono] ^cosso cpuitt 

43frtS^» j <Dm& "7" tcJlLTt^5b^5^«:^x 
7^t§o ftt, *^»j<Dfit# "7" tc3SLTl^ 

ftv^s^tcti. cpunt s 5 o 8izi3^rm& j * 
[0111] ctitcMLTs mm£jzmfz>^mj <om 

ifi "7" icmhXh^tS 5 0 7 iCTf mLfcmiklc 
fin C P U 1 #aJI£: S 5 0 9 tc)ia6£ 0 CO S 5 0 so 
9 CPUUfc S«FjSk:fefrt«^R i ©fflKA 

"7" ICiLTt^*^5fr*fi^*no ?LT, 
*r£^»i<Dffi;W "7" tcllt^^^i'&^ti, CP 
U Hi, SSlOtio^Si^^iJ^^Uft 

1SS 5 0 5 -\M*To 

[0 112] £*UC*fLT, ffl^tc*5t3-§^» i Offi 
W "7" fcjSLTV^^S 5 0 9lcm^Lfci§^lz 
fcJ\ CPU Hi, £QJI£*S 5 1 l'NitfeSo CCDS 5 1 

i -ea, c p u ni, s 5 o 6 tcrwrn^nfc ixj f 40 
w%m» mi 0 (d) ) tiwi-©iB«T f -^ (0 

11 (e) ) £&£ 0 



[0 113] 

H#*P 5 ttUSO^ nlflga^SlT?, JSBflfrr- * t£ f$gijfif 
£0 ti£oT. B«7^#<DWflJ#*fcJ\ HCBfc-r-** 

• gar Si^s^av^T, ASMcoi2iis«^ffl 
[0 1 ] ^mmmmm 

[02] ©*tt»j»tccfe«siiffl3>tf 

i - # £>«WM»£«:*rr :/ n ^ 0 
[0 3] 0 2OC PU 1 KTHff*hSllBiJW«a3a 

[04] *mmm 1 (DmMmmic&%M&m^^¥ 
[05] 04ocpui icrmni£H2>mmmmt&& 
[06] mmm 1 tc<fc«»But»«aa«Bi<DrtSF*^ 

[0 7] 06CDS2O 4tCT^T*nS2^7£DCT 
[0 8] 06 OS 2 1 3tCTHtf^n5 2^7Ca?DC 

[0 9 ] mmm 1 k ± s«guffifflaftm«ia<Drts*^ 

[0 l 0 ] ^JfiM l lc £ £*£B^i£ffiS<DSKn«: 

[0 1 i ] *fifc#!i i £ jcss£B«wwffl3^ft* 
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SB 
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4 1- 
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[H3] 

\2<dc pu 1 izTnftznzm%mmm&&&cDft&&^T7n--zhv- h 



S 0 1 0 



r=i + 1 



S 0 0 1 



S 0 0 2 



S 0 0 3 



S 0 0 4 



S 0 0 5 



s o o e 



s o 0 7 



fjCcj^Sj^iUfc-riq m <D 
«C lk (l^KSKj)«r*^ 

3 



S O 0 8 



N O 




S0O9 



soil 



SOI 2 



^ END ^ 
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[04] 



2 2- 

3 4- 
3 3- 



mm 



ROM 



CPU 



1 6 



1 7 



1 8 



C87] 



S 3 1 1 



W i (X)= V n-«* \m 



:^7TDCT 



S 3 0 1 



i =0 

3EI 



S 3 0 2 



S 3 0 3 



i =0 



S3 0.4 



1 S 3 08 ^ s 305 






C(ij)= "z z'GCx,y)-w,(x)-Wj(y) 

x-O y-O 




s 3 0 7 — -^-^ S 3 0 fi 


j=i+i 





S 3 1 0 



.YES 



N O 



S 3 0 9 



(17) 



# 10-234012 



[0 5] 



S 1 0 8 



im i o] 



CEfEED 



S 1 0 2 



S 1 0 3 



S 1 0 4 



S 1 0 5 



I = 1 



S 1 0 6 



Sj -f, (C,) 



N O 




1 0 7 



S 1 0 9 



END 



(a) 

MM 



2ftxDC7 




[3 



O 8 ) t*HPF 1 0-2 3 4 0 1 2 



CH6] 



S 2 0 1 



S 2 0 2 



S 2 0 3 



[x^-(n-i),v=o^<n-i)] 



S 2 0 4 



J£ffltt£ 2^7tDC T&fcr £ 



S 2 0 5 



f =0 



S206 



j = 0 

5 



S 2 0 7 



*c iJk )-S<i,j)**fc«ija(D 
*¥C ijk (1£K£Kjj)*#«>5 



S 2 0 8 



C(i , j)£Cj jk **t>C(i J)fCi5L\J» 



j - j + l 



S 2 0 9 



S 2 1 2 




S 2 T 1 



YES S 2 1 3 



MStt*HliMC(i,j) [MMNM>, 
j=0-CN-l)] £2;*5ti£DCT 



S 2 1 4 



*«SMi&*MR()(,y) 

[x-CK-(N-I).y=0-(l^l)J 

^ END ^ 



WM¥ 10-234012 



[H8] 



x-x+1 



S 4 0 1 



Wo(x)=/L i Wi(x)=V^-cm "< 2 2 ~ T i 



i. 



S 4 0 2 



x = 0 



S 4 0 3 



y = o 



S 4 0 4 



R(x f y)=Z |c(i,j)-Wi(x).Wj(y) 



S 4 0 5 




Mil] 




(20) 



t#HS¥ 1 0-2 3 4 0 1 2 



[139] 



S 5 0 1 



S 5 0 2 



»»HMHfcR(x. y) 

tx=0~(N-l)/M)MN-1)] 



S 5 0 3 



S 5 0 4 



i = 0 



5 



S 5 0 5 



j =0 



S 5 0 6 



s r (ij)=f(c(ij))*jtt±i-r-£ 



S 5 0 8 



j = J + 1 



N O 



S 5 0 7 



S 5 1 0 



= i+1 



NO 




S 5 0 9 



S 5 1 1 



S'(iJ) [i-0-7, j=(W] 



END 
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# HBT 1 0-2 3 4 0 1 2 



(54) mm^—^^nmrnmnm^^m, mfflmmm®&%tircwm7 s ~-zfrz<Dm,wmnmBi 



